30th International Cosmic Ray Conference 




ICRC*/ 



Merida, Mexico 



o 
o 

(N 



o 



Oh 
i 

o 



> 
in 

o 
o 



X 



Relative abundances of cosmic ray nuclei B-C-N-O in the energy region from 10 
GeV/n to 300 GeV/n. Results from ATIC-2 (the science flight of ATIC). 
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Abstract: The ATIC balloon-borne experiment measures the energy spectra of elements from H to Fe 
in primary cosmic rays from about 100 GeV to 100 TeV. ATIC is comprised of a fully active bismuth 
germanate calorimeter, a carbon target with embedded scintillator hodoscopes, and a silicon matrix that 
is used as the main charge detector. The silicon matrix produces good charge resolution for protons and 
helium but only partial resolution for heavier nuclei. In the present paper, the charge resolution of ATIC 
was improved and backgrounds were reduced in the region from Be to Si by using the upper layer of 
the scintillator hodoscope as an additional charge detector. The flux ratios of nuclei B/C, C/O, N/O in 
the energy region from about 10 GeV/nucleon to 300 GeV/nucleon obtained from this high-resolution, 
high-quality charge spectra are presented, and compared with existing theoretical predictions. 



Introduction 

The ATIC spectrometer, its calibration and the al- 
gorithm of trajectory reconstruction have been de- 
scribed [3, 9, 7]. Charge resolution provided by 
the silicon matrix is sufficient to obtain spectra of 
primary protons and helium [8, 6] and preliminary 
spectra of some abundant heavy nuclei [5, 6]. 

Very important to understand the mechanism of 
propagation of cosmic rays in the Galaxy is the 
boron (which is a secondary nuclide) to carbon 
ratio in cosmic rays. The problem of B/C ra- 
tio has been experimentally investigated in the en- 
ergy range 0.5-50 GeV/n (see [1] and references 
herein). The energy range of the ATIC experi- 
ment allow data for higher energies (up to 200- 
300 GeV/n) to be obtained. But there are obsta- 
cles: 1) low charge resolution of the silicon matrix 
in the range of charges 5-8 and 2) high background 



in the silicon matrix charge spectrum in the range 
of boron and carbon (see fig. 1). In this paper we 
use the upper layer of the scintillator hodoscope to 
improve the charge resolution and to reduce back- 
grounds in B-C region to measure B/C in the ATIC 
experiment. 

Improved charge spectrum 

The upper scintillator layer of the hodoscope is 
comprised by 42 parallel scintillator strips 1 x 2 x 
88.2 cm 3 . Using these scintillators as a supplemen- 
tary charge detector, involves a multi-step proce- 
dure of calibration and normalization of the sig- 
nals which will be described in detail elsewhere. 
In brief, the method is the following. The first step 
is to use the usual method to measure the charge of 
primary particles - trajectory reconstruction from 
signals in BGO-calorimeter project - to the silicon 
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Figure 1 : The charge spectra obtained with the sili- 
con matrix only and with the silicon matrix plus the 
upper layer of hodoscope for the range of energy 
deposit in BGO calorimeter 50-100 GeV (primary 
energy per particle approximately 150-300 GeV) 



matrix. The charge detected in the silicon matrix 
(QsO is the maximal signal in the area of confu- 
sion for the trajectory [3, 9, 7]. In the second step, 
we find the charge detected in each scintillator strip 
and select the strip with the charge nearest to the 
charge detected by the silicon matrix (Qsd)- This 
charge is accepted if the distance from the strip to 
the reconstructed trajectory is less than 5cm and 
is rejected otherwise. In the third step, an event 
is rejected if |Q Sci — Q si \ > 0.25. The final re- 
sult for the charge is Q = (Qsd + Qsi)/2. This 
procedure reduces the backgrounds in the charge 
spectrum and improves its resolution, but reduces 
the initial statistics by a factor of about 4. There 
are other strategies to process charge data from the 
scintillator hodoscope, but for this paper we se- 
lect the strategy of "higher resolution - lower back- 
grounds - lower statistics". The charge spectra ob- 
tained with the silicon matrix only and with the sil- 



icon matrix plus the upper layer of hodoscope are 
compared in fig. 1. 



Measurement of the relative fluxes 

We calculate the ratio of fluxes of different nu- 
clei in cosmic rays to the flux of carbon against 
energy of particles per nucleon. It is a multi-step 
procedure which is designed to obtain the most ex- 
act information for fluxes of nuclei with charges 
4 < q < 14. 

1. For five ranges of the energy deposit Ed in 
the BGO calorimeter (50-100, 100-200, 200-398, 
398-794, 794-1585 GeV) we obtain the charge 
spectra (similar to fig. 1, lower graph), and decom- 
pose each by Gaussian fits (the value \ 2 per degree 
of freedom is close to 1 in all cases). The positions 
of peaks are determined and charge cuts are devel- 
oped for each particular primary particle such that 
the margins of cuts are at the half of path between 
adjacent peaks. The number of counts I® q in each 
charge bin q for Ed range number s is the raw data 
to obtain the fluxes of primary particles (s = cor- 
responds to the energy region of Ed 50-100 GeV, 
etc). 

2. Protons and helium interacting in the material 
(aluminum honeycomb and other) of ATIC above 
the silicon matrix can sometimes simulate heav- 
ier nuclei. This effect is energy dependent (grows 
with energy). Corresponding backgrounds Bs'^ e 
for each value 7, are calculated by simulation 
of propagation of protons and helium through the 
ATIC instrument by the FLUKA code [2], with 
simulation of the conditions of charge selection 
(see previous section). The apparatus charge line 
widths are accounted for as well. This proce- 
dure produces the corrected values of intensities 
l] „ — 7?„ — Bs'^ e . The value of Bs'^ e for boron 
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(q = 5) varies from 9% to 36% of I° s q . 

3. Particles with charges q > 15 fragmenting in 
the material above the silicon matrix can also pro- 
duce nuclei of 4 < q < 14. The corresponding 
backgrounds were subtracted but the effect is small 
(about 0.1% for boron) and we do not describe the 
method of subtraction here. 

4. Each nuclei of 4 < q < 14 due to interac- 
tions in ATIC, and due to the apparatus broadening 
of the peaks, produces a "charge response" of the 
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device which may be described for each particu- 
lar nucleus q at the entrance to the instrument by a 
set of the coefficients K\ , K q , . . . , Kf 4 , where K\ 
is the probability to find nucleus q in the charge 
bin 4, etc. Of course the coefficient dominates 
strongly in the set K\, . . . , Kf 4 . In other words 
the matrix | \K q , \ | is diagonally-dominated. Let F q 
be the intensity of the nucleus q at the entrance to 
ATIC. Then, for each energy region s, the experi- 
mental charge spectrum J* after subtraction of p- 
He backgrounds (described above; here, we do not 
write the index s for simplicity) may be written as 

I\ = K\F 4 + KlF 5 + ... + K\ A F 14 

(1) 

l\ 4 = K\ 4 F 4 + K* 4 F 5 + ... + K\\F 14 

Eq. (1) is a square linear system relative to un- 
known values i<4, F 5 , . . . , F14 with diagonally- 
dominated matrix elements, and it can be eas- 
ily solved by usual methods. The coefficients of 
the system K q , are calculated by simulation of 
propagation of different nuclei through ATIC with 
FLUKA and the values I q are already known af- 
ter steps 2 and 3. The ratio of B/C (calculated as 
the ratio of the contents of the related charge bins) 
reduced by \\%-A2% for different energies. 

5. The next step is a transition from the spectra of 
fluxes as a functions of Ed (the result of step 4) 
to the spectra in primary energy per nucleon. For 
each nucleus this procedure includes firstly a cal- 
culation of expected primary energy for each edge 
of the region of Ed (see step 1). To solve this prob- 
lem FLUKA simulation of energy deposition was 
used with supposition of the primary differential 
momentum spectra to be a power law with index 
7 = —2.6 (there is only a weak dependence of 
exact value of 7). After normalization of the pri- 
mary energy to the atomic weight and to the width 
of the related energy region, the fluxes of different 
nuclei have different energy per nucleon binning. 
To obtain the ratio of fluxes at the same energy we 
calculate the energy points obtained as a geomet- 
rical mean for the corresponding points for boron 
and carbon and calculate all fluxes for these energy 
points by interpolation of the spectrum of each nu- 
cleus. This procedure increases B/C ratio from step 
4 by 13%— 23% (different for different energies). 

6. The mean altitude of the flight of ATIC-2 was 
36.5 km which corresponds to 4.87 g/cm 2 of resid- 



ual atmosphere. To obtain the primary fluxes above 
the atmosphere, the interaction of nuclei in the at- 
mosphere should be accounted for. The interaction 
may be described as fragmentation of nuclei with- 
out changing the energy per nucleon. Then the in- 
teraction for each primary energy and for each pri- 
mary nucleus q may be described by a set of co- 
efficients L q q , ,q'<q which show the probability 
to find the nucleus q' at the entrance of the instru- 
ment. The coefficients L q , were calculated by sim- 
ulation of propagation of nuclei in the atmosphere 
by FLUKA. Let ip q be the flux of nucleus q in en- 
ergy per nucleon at the entrance of ATIC (these 
values are known after step 5; we omit the index 
of the energy for simplicity) and tp q be the same 
values above the atmosphere for some definite en- 
ergy per nucleon. Then one can write 

V>4 = Ljip 4 + ... + L\ 4 (p 14 + e 4 
(2) 

1pl4 = L\\Lp 14 +£14, 

where e 4 , . . . , £14 are small corrections related to 
the fragmentation of nuclei heavier than silicon 
(q = 14). If s q are known then the system (2) is 
a square linear system relative to ip4, . . . , ipu with 
the triangle and diagonally-dominated matrix L q ql 
and it can be solved directly starting from the final 
equation. The atmospheric correction reduces B/C 
ratio by 13%— 33% for different primary energies. 

The experimental errors were calculated as a com- 
bination of the Poisson dispersion of the exper- 
imental statistics and the statistical errors of the 
simulations. If the desired quantities were obtained 
as a solution of a linear system (as in steps 4 and 
6) then corresponding complete covariation matrix 
were calculated by the Monte Carlo method. All 
reported errors are standard deviations. 

Results and discussion 

The results for B/C, N/O, C/O ratio are presented 
in table 1. The data for B/C, N/O and C/O ra- 
tio along with the data of HEAO-3-C2 experiment 
[1] with theoretical predictions are shown in fig. 2. 
One can see that the data of present work for B/C 
and N/O is somewhat above the data of [1] but is 
extended to higher energies. The theoretical curves 
in fig. 2 are calculations in leaky box approxima- 
tion. The dashed line is based on the HEAO-3-C2 
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Table 1: B/C, N/O, C/O ratios as a function of pri- 
mary energy (GeV/n). The numbers in parenthe- 
sis give the uncertainty in the last significant digits 
quoted. 



E 


B/C 


N/O 


C/O 


19.9 


0.180(11) 


0.219(10) 


1.020(26) 


38.3 


0.169(15) 


0.199(13) 


1.087(43) 


74.3 


0.119(29) 


0.184(24) 


0.933(60) 


149 


0.156(53) 


0.172(39) 


0.934(105) 


307 


0.064(63) 


0.144(68) 


1.022(227) 
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Osborn & Ptuskin, leaky box model [4] 




HEAO-3 model, leaky box model [1] 
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fit for the Galaxy escape length [1] A esc = 34.1/3i?-°- 60 
(R is the rigidity) and the solid line is for the escape 
length obtained in the model of Kolmogorov type 
of magnetic turbulence and reacceleration during 
propagation [4]: 

Aesc = 4.2(i?/i? )- 1/3 [1 + (R/R a )- 2 / 3 ] gem- 2 , 
where i?o = 5.5 GV. Whereas the experimental 
data support general trend of decreasing B/C and 
N/O ratio with energy, it is impossible to distin- 
guish between different models of propagation of 
particles due to the experimental uncertainties. 

It should be noted that our experimental data are 
model dependent due to extensive simulation of 
the backgrounds by the FLUKA code, but this 
could be improved by usage of additional simu- 
lation codes. One can expect that the experimen- 
tal separation between different models of propa- 
gation would be possible with additional experi- 
ments. 
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Figure 2: B/C, N/O and C/O (top-down) ratio from 
this work, from HEAO-3-C2 experiment and from 
leaky box models. 
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